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Abstract 

In this note, we solve the Loewner equation in the upper half- 
plane with forcing function for the cases in which £(t) has a 
power-law dependence on time with powers 0, 1/2 and 1. In the first 
case the trace of singularities is a line perpendicular to the real axis. 
In the second case the trace of singularities can do three things. If 
£(t) = 2\/~Kt, the trace is a straight line set at an angle to the real axis. 
If = 2a/ k(1 — t), as pointed out by Marshall and Rohde the 
behavior of the trace as t approaches 1 depends on the coefficient k. 
Our calculations give an explicit solution in which for k < 4 the trace 
spirals into a point in the upper half-plane, while for k > 4 it intersects 
the real axis. We also show that for k = 9/2 the trace becomes 
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a half-circle. The third case with forcing = t gives a trace that 
moves outward to infinity, but stays within fixed distance from the real 
axis. We also solve explicitly a more general version of the evolution 
equation, in which £(t) is a superposition of the values dbl. 

Key words: conformal map, Loewner equation, singularity, slit 
domain. 

1 Introduction 

The Loewner differential equation was introduced by Karl Lowner ^T] (who 
later changed his name to Charles Loewner) to study properties of univalent 
functions on the unit disk. The differential equation is driven by a function 
that encodes in an ingeneous way continuous curves slitting the disk from 
the boundary. A few years ago, Oded Schramm realized that by taking 
Brownian motion as the driving function, the Loewner equation generates 
families of random curves with a conformally invariant measure, that he 
called "Stochastic Loewner Evolutions" (SLE's) [THj . Moreover, he showed 
that if the loop-erased random walk has a scaling limit, and if this scaling 
limit is conformally invariant, then it must be described by an SLE, and he 
made similar conjectures about the scaling limits of uniform spanning trees 
and of critical percolation. The existence of conformally invariant scaling 
limits and the connections with SLE were established for the first two models 
by Lawler, Schramm and Werner j^j, and for critical site percolation on 
the triangular lattice by Smirnov |14j . Over the past few years Schramm's 
idea has been investigated further, mainly in a combined effort of Lawler, 
Schramm and Werner jS], leading to new developments in the study of (the 
scaling limit of) discrete models. 

In SLE an ensemble of driving functions generates an ensemble of shapes. 
It is a difficult task to understand fully the relation between these two ensem- 
bles. For this reason we investigate here the simpler problem of the relation 
between a driving term in Loewner's equation and the curve it generates, by 
explicit calculation. In particular we study exact, deterministic solutions of 
Loewner's equation in the upper half 7i of the complex plane. One of the 
purposes of this study is to present some explicit examples that may help to 
elucidate the behavior of the curves generated by Loewner's equation. For 
example, our solution for square-root forcing in section 0] shows a transition 
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Figure 1: This figure shows how a slit 7 in the z-plane maps into a linear 
boundary in the w-plane. Notice that two neighbouring points in the z-plane 
have images which are far apart in the w-plane. 



from a simple curve that stays above the real line to a curve intersecting the 
real line, which is similar to the transition occurring at k = 4 in SLE. 

To get us started, we need to understand how the Loewner equation 
describes curves slitting the half-plane Ti. So suppose that j(t), t > 0, 
is a simple curve in TC emanating from the real line. Then according to 
conformal mapping theory, there exists a unique conformal mapping of the 
form w = gt{z) that takes the domain 1Z consisting of all points in 7i minus 
those on the slit 7[0, t] onto the upper half of the w-plane, in such a way that 
near infinity this map has the expansion 

g t (z) = z + c(t)/z + 0(l/z 2 ). (1) 

See figure ^ for a schematic picture of the mapping. It is a theorem that 
the coefficient c(t) (which is called the capacity) is continuously increasing 
with t. Therefore, the parameterization of the curve can be chosen such that 
c(t) = 2t, and henceforth we assume that this is the case. 

The conformal maps g t satisfy a surprisingly simple differential equation, 
which is the Loewner equation 

dt~ gt - m { } 

with the initial condition 

go = z (2b) 

for all z in Ti. Here, the value of £ at time t is just the image of the point 
7(t) under the map g%. We call the driving function or forcing function 
of the Loewner equation. It is continuous and real-valued. (For a text-book 
discussion of Loewner's equation and proofs of the above statements when 
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the domain is the unit disk, see [2] and [3]. The half-plane case we are 
considering here is analogous and is discussed in [7] and jS]). 

Conversely, the Loewner equation (j2J) has a solution for any given contin- 
uous real- valued function £(;£). This generates the function gt(z). For each 
value of t this function can be thought of as a mapping of the form w = gt(z) 
which takes some connected subset of the upper-half z-plane, 71, into the 
entire region 7i above the real axis of the w-plane. Correspondingly, there 
exists an inverse function ft(w), which obeys g t (ft(w)) = w for all w in the 
upper half-plane. This function maps the region 7i of the w-plane into the 
region 1Z of the z-plane. It obeys the partial differential equation 

df t (w) _ 2 df t (w) 
dt w — £ (t) dw 

The Loewner equation continually generates new singular points of gt 
that are mapped by g t onto the corresponding points £(£) in the w-plane. 
Thus, as time goes on, points are continually removed from TZ. If £(t) is 
sufficiently smooth [T2], these singular points form a simple curve z c (t), where 
the points z c {t) obey 

gtMt))=t{t). (4) 

This curve is called the trace of the Loewner evolution, but we also refer to 
it as the line of singularities. (Naturally, if we take to be the driving 
function corresponding to a given curve j(t) as described above, then z c (t) 
will coincide with 

Equations (J2J) have a simple scale invariance property under changing 
both the scale of space and of time. Any change of time of the form 1 1— > a 2 t 
can be compensated by a change of scale. To do this, define a new g and £ 
by 

9t(z) ■= -g a 2 t (az) and = -£(ct 2 t). (5) 

a a 

This pair then satisfies equations (J2J). Note that by the Schwarz reflection 
principle of complex analysis pQ, the map g t extends to the lower half-plane 
and satisfies gt(z) = gt(z). It follows from equation (jSJ) that flipping the sign 
of the forcing has the effect of reflecting the regions 1Z and Ti, (and hence 
also the line of singularities z c (t)) in the imaginary axis. Likewise, a shift in 
the driving function £(t) can be compensated by a similar shift in g t , since 
the pair 

g t (z) := g t (z - a) + a and = + a (6) 
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again solves equations (J2J). This shows that a shift in the driving function 
simply produces a shift of the same magnitude in the trace z c (t). 

In this note we are interested in the question how the behavior of the trace 
is related to that of the driving function. This question can be investigated 
both numerically and, in some cases, by exact methods. A simple numerical 
method to find the trace generated by a given driving function is to nu- 
merically integrate Loewner's equation backwards from the initial condition 
9t = to obtain the value of go = z c (t). To deal with the singularity in 
Loewner's equation, we can use the first-order approximation 

2 dt 

9t - gt-dt H T77 T7v {<) 

9t-dt -£(t- dt) 

Substituting g t = £(t) and solving this equation for gt-dt gives an approxima- 
tion for the value of g at time t — dt, from where we then integrate backwards 
to obtain g . This method is good enough for the purpose of producing pic- 
tures. The figures El IU El and [7| in this note were obtained in this way. 

Exact solutions of equations (J2J) can be found in a few cases. The simplest 
are ones in which the forcing has one of the forms 

i{t) = C{l-tf (8a) 



or 

£{t) = Ct 13 (8b) 

for (3 = 0, 1/2, and 1. For a driving function of the form of equation (j8b|) . a 
rescaling according to (JSJ) gives the new forcing 

= Co? p ~H p . (9) 

This result allows us to scale away the constant C when (3 ^ 1/2. For the 
special value (3 = 1/2 the multiplicative constant can not be scaled away, 
and can be important in determining the form of the solution. 

In sections El El and 0] we shall derive the solutions respectively for the 
powers 0, 1, and 1/2. Section El contains an explicit construction of the 
maps g t for a half-plane slit by an arc, and we point out that this is a special 
case of square-root forcing with a finite-time singularity. In section El we 
consider a more general version of Loewner's equation, where the forcing is 
a superposition of the values ±1 and produces two lines of singularities. 



5 



A 



wq- A w + 



Figure 2: The trace produced by a constant forcing is a vertical line segment. 

2 Constant Forcing 

Here we look at the almost trivial situation in which {3 = and we can take 

H(t) = A. (10) 

This situation is well-known (see for example |2| and appendix A in jH]), but 
for the sake of completeness we also treat it here. We could of course shift A 
away by using equation ©, but this case is so easy that shifting is hardly 
worth the effort. The equation for g can be solved by inspection, giving 

g t {z) = A+[{z-Af + At] 1 /\ (11) 

The inverse map is of the same form as the direct one, namely 

f t (w) = A+[(w-A) 2 -At} 1 / 2 . (12) 

At time t, the function g acquires a new singularity at the point 

Zc (t) = A + 2it 1/2 (13) 

which maps under g t to the point £(t) = A in the u>-plane. Thus the singu- 
larities form a line segment parallel to the imaginary axis and the mapped 
region of the z-plane is the upper half-plane minus that line segment. The 
mapping has a singular point at the place where the trace meets the real axis. 
We can approach this point z = A from two different directions, giving us 
two different images at the points 

w ± = A±2t 1/2 . (14) 

Figure El gives an illustration of the mapping. 

Because any smooth enough function £(t) looks locally as if it were con- 
stant, we might expect that smooth £(t)'s generate traces of singularities 
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which look like simple curves coming up from the real axis. In fact, one can 
prove that this statement is correct whenever is sufficiently smooth so 
that it is Holder continuous with exponent 1/2 and sufficiently small norm, 
i.e. so that 

^- S) - m <C (15, 

for all t and some constant C. The tricky issue is, when can it be that the 
curve will touch itself or hit the real line. Condition (j!5|) is enough to ensure 
that this will not happen, as was shown by Marshall and Rohde [12] , who 
did not give the value of the constant C. Our solution in section 0] below 
gives a specific example of the transition occuring at C — 4. In a parallel and 
independent paper, Lind ^U] proves also generally that H61der-l/2 functions 
with coefficients C < 4 generate slits. 



3 Linear Forcing 

The next case, not much harder, has the forcing increase linearly with time 
according to 

m = t. (i6) 

We do not need an additive or multiplicative constant in equation (|16|) be- 
cause these can be eliminated by equations (j5J) and (JEJ). If we now redefine 
the independent variable to be h = g — £, then h obeys 



dh 2 — h dh 



(17) 



dt h dF{h) 
where F(h) is the function 

F(h) = h + 2\n(2- h). (18) 
In terms of F the solution is 

F(h) = -t + c(z). (19) 

Here, c(z) is a constant which must be set from the initial condition, g = z, 
giving the solution F(h) = — t + F(z) or equivalently 

F(g-t) = F(z)-t. (20) 
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Figure 3: The shape of the line of singularities in the z-plane that arises from 
the forcing £(£) = t, up to time t = 100. 

To find the trajectory of the singularity, use equations (jlj) and (|2U|) to get 

F(z c {t)) = F(0)+t = 2m2 + t (21) 

It is possible to find an explicit expression for this trajectory. To do so, 
substitute 2 — z c (t) = r t exp(—i(j) t ) into the previous equation. Then split the 
equation in its real and imaginary parts to obtain 

21nr t -r t cos(j) t = 21n2 + t - 2, (22a) 
r t = 2&/sin&. (22b) 

After substituting the second equation into the first, it can be shown that </> 
increases monotonously in time from the value 0o = to (f)^ = it. 

In terms of the parameter 0, the line of singularities may be written 
explicitly as 

z c (t) = 2 - 20t cot fa + 2i(j) t . (23) 

This shows that the line of singularities moves outward to infinity while 
remaining within a fixed distance of the real axis. Figure El shows the first 
part of this trajectory. 

For small and large values of t we can also construct asymptotic forms of 
the solution (|2*T|). As t goes to zero, one can expand F(z c (t)) to third order 
in z c (t) and find that, to this order, the result is 

z c (t) =2it 1 ' 2 + -t + 0(t z/2 ). (24) 
3 
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For large t, 

z c (t) =t-21n[(t-2)/2] + 2vri + 0(lnt/t). (25) 



4 Square-root Forcing 

The next case, much more interesting, has the forcing be a square-root func- 
tion of time. This case has two subcases: one in which the forcing has a 
finite-time singularity based on the rule of equation (J8a|) . i.e. 



£(t) = 2[k(1 -t)] 1/2 , t<l, k>0. (26) 

Here, as we shall see, the result changes qualitatively as k varies, with the 
critical case being k = 4 where the singular curve intersects the real axis. In 
the other situation, the forcing has an infinite-time singularity with the rule 

£(t) = 2[Kt} 1/2 , k>0. (27) 

In this case z c (t) is just a straight line, as is already clear from the scaling 
relation ®. 

4.1 Infinite-time Singularity 

To find the solution for the forcing (|27jh define the new variable G = g/t 1 ^ 2 
and set r = In t. Then G satisfies 

^ = - G/2 + 2/(G -^ ) = (^-»-) <*) 



where j/ ± = k 1/2 ± (k + 4) 1/2 . It follows that 

y+ 



dG 
dr 



G-y + G -y_ 



\{y+-y-)- (29) 



Therefore, if we set 



;= 2»>MC-,-)-2,-h(C-, t ) 

then dH(G)/dr = — 1 which integrates to 

-if (GO =r + c(z) = In t + c(z). (31) 
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The constant c(z) can be determined from the observation that in the 
limit when t approaches 0, 

H(G) +\nt = 2 y+ ln (9-y -t 1/2 ) - 2y- Info ~ V+t 1/2 ) 



y+-y- 

— ► 2\nz. (32) 

Our solution (j31)) for general t becomes simply 

H{g/t l / 2 ) = 2\n{z/t 1 / 2 ). (33) 

Since the line of singularities is determined by the condition that g is equal 
to the forcing we get 

z c (t) = Bt 1/2 where B = exp[-H(2K 1/2 )]. (34) 

More explicitly the coefficient is 

1 K V2 



B " 2( (« + 4)V2-«i/2j 6XP 



1 Yi 

7TZ 1 



(35) 



2 V (k + 4) 1 /2 
so that the line of singularities is set at an angle to the real axis which is 

1 / K 1 / 2 \ 

= 7F 1-7-^ • (36) 



2 V (^ + 4) 1 /2 

For k = the line is (as we know) perpendicular to the real axis while as 
n — > oo the angle of intersection becomes smaller and smaller. 

4.2 Finite-time Singularity 



Now we turn to the forcing (|2o|) with a singularity after a finite time. In this 
case, one can eliminate the time- dependence by using the new variable 

G = g/{l-tfl 2 . (37) 

Then G obeys 

An 

^ = G/2 + 2/{G-2k 1 ' 2 ) (38) 
dr 
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where r = — ln(l — t). Once again the derivative becomes a simple function 
of the unknown, again appearing in the form of a ratio of polynomials, i.e. 



_dG (G-y + )(G-y-) 
dr 2(2/^/2 _G) • ^ 

where this time the roots are y± = n 1 / 2 ± (k — A) 1 / 2 . Notice that now the 
roots are real and positive when k > 4 but that they are complex for k < 4. 
Thus we can expect a qualitative change in the solution at k — 4. 

The integration of the equation proceeds exactly as in the previous case, 
giving the solution 

H(G) = t + c(z) (40) 
where H has the same form as before, but with different y±: 

H(G) = 2 y+HG-y-)-2y_HG-y + ) ^ 

y+-y- 

The initial condition then implies that 

H( S /(l-t) 1 ' 2 ) = -Hl-t)+H(z) (42) 

and the equation for the line of singularities is then 

H(2k 1/2 ) = -ln(l -t) + H(z c (t)). (43) 

Note that equation (|4Hjl for the line of singularities tells us that H(z c (t)) must 
approach — oo when t approaches 1. In the following subsections we shall 
consider the asymptotics of z c (t) in this limit in more detail. We separate 
the discussion into the two cases « < 4 and k > 4. Then the section is 
completed with the derivation of the shape of the critical line of singularities 
at k — 4. 
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4.2.1 The logarithmic Spiral for k < 4 

When k < 4, the real and imaginary parts of H(z) are 

k 1/2 

ReH(z) = ^ - /2 Arg (z - y_) + In \z - y_| 



~ (4 _ K )i/2 Arg ~ y+) + ln l z ~ ( 44a ) 

Imif(z) = Arg (z - y_) - -—- 2 ln \z - y_| 

« 1/2 

+ Arg (z - y + ) + _ /2 ln |z - y + |, (44b) 

from which it follows easily that H(2k 1 ^ 2 ) is a positive real constant. Since 
the real part of H(z c (t)) must go to — oo when t approaches 1, we conclude 
that z c (l) = y + , since the other possible candidate for the limit, i.e. y_, is in 
the wrong half-plane. It then follows from the observation that Im H(z c (t)) 
must vanish, that the trace must have wrapped around y + an infinite number 
of times. Hence the line of singularities is spiralling in towards the point y + , 
as is depicted in figure |U 

To find a more explicit expression for the asymptotics, notice that when 
z c (t) is close to y + , we can make the approximations 

\n\z c (t) -y-\ ta ln (2(4- k) 1/2 ) and Arg (z c (t) - y_) w vr/2. (45) 



Splitting equation (}45|) in its real and imaginary parts, and substituting the 
above approximate values then gives a system of two equations that we can 
solve for the unknows ln|z c (£) — y + \ and Arg (z c (t) — y + ). This gives the 
result 

\n\z c {t)-y + \ w [A(«) + (4-k)1h(1-*)]/4, (46a) 
Arg(z c (t)- y+ ) « [ J B( K )- K 1 / 2 (4-K) 1 / 2 ln(l-t)]/4, (46b) 

where the constants are 

A(k) = In 16 + (k - 2) ln(4 - «) 

+ (2Arg(y + )-7r)/t 1/2 (4-/t) 1/2 , (46c) 
= K 1 / 2 (4-K) 1/2 ln(4-K) 

+ (tt - 2 Arg(y+))K - 2tt. (46d) 
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-2.0 -1.5 -1.0 -0.5 

Figure 4: The trace of singularities for square-root forcing with different 
values of k < 4. Also shown is the limiting case k = 4 where the trace starts 
to hit the real line. All the curves have been translated to make them start 
in the origin at time t = 0. 

Thus, the distance between z c (t) and y + is decreasing like a power law in t, 
whereas the winding number of z c (t) around y + grows only logarithmically. 

4.2.2 The Intersection for k > 4 

For k > 4 the real and imaginary parts of H(z) are 

+ V 1 " (k- 4)V2 jlnk-y+l, (47a) 
Im#(z) = (l+ (/t ^ 2 )1/2 )Arg(^-y-) 

+ ^ ~ (K _ 4) i/ 2 jArg(s-y + ), (47b) 

and H(2k 1 ^ 2 ) is again a positive real number. This time, as t approaches 
one, z c (i) must approach the point y_ on the real line, see figure El In fact, 
we can calculate the angle at which the line of singularities intersects the real 
line from the observation that Im H(z c (t)) must vanish. 



Re H(z) 
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Figure 5: For a square-root forcing with k > 4 the line of singularities 
intersects the real line at an angle that varies with k. The curves were 
translated to make them start in the origin. 



Indeed, if we denote by <fi the angle at which z c (t) hits the real line 
(measured from the real line in the positive direction), then we have that in 
the limit t — > 1, Arg(z c (t) — y + ) —* it whereas Arg(z c (t) — yJ) — > <p. The 
condition Im H(z c (t)) = then gives 



1/2 _ ( K _ A \l/2 



71 



4) 1 



«V2 + 



K 



4)V2 



(48) 



when t goes to 1. As we can see, we have a glancing incidence for k = 4 and 
K — ► oo, while for k = 9/2 the incidence is perpendicular. In fact we can 
prove that for k = 9/2 the line of singularities is a half-circle of radius y/2. 
We will do this in section where we give an explicit construction of the 
conformal maps g t when the half-plane is slit by a half-circle. 



4.2.3 The critical Curve at k = 4 

It is not difficult to find an explicit expression for the critical line of singu- 
larities at k — 4. In this case, the two roots in equation (|39|) are the same, 
y± = 2, and the equation can be rewritten as 



dG 



G-2 (G-2) 2 J dr 
Integration is straightforward and gives 

4 



1. 



21n(G-2) + 



G 



ln(l-t)+21n(z-2) 



(49) 



(50) 
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Figure 6: The construction of the map ip s taking Ti minus the slit 7[0, s] 
onto H, in three steps. The illustration also tracks the images of the upper 
half of the unit disk and of the four "special" points —1, 0, j(s) and oo. The 
definitions of j(s), R s and a s are in the text. 



where we have determined the integration constant from the initial condition 
Co = go — z. The line of singularities is determined by the condition that 
for z = z c (t), G = 4. 

Substituting z c (t) — 2 = r t exp(z0 t ) and splitting the equation in real and 
imaginary parts leads to the expression 



sin 2cf) t . sin 



2 



z c (t) = 2 + — + 2i— (51) 

<Pt <Pt 

for the critical curve in terms of the parameter <fi. This parameter <ft increases 
monotonously with time from </> = to <fii = it. 



5 Loewner Evolution for a growing Arc 

In this section, we consider the situation where the upper half-plane TC is slit 
by the curve 7(5) = exp(is), s G [0, 7r). The idea is to construct the conformal 
maps that map the half-plane minus the slit 7(0, s] back to the half-plane 
and have the expansion of equation (JIJ). After reparameterizing time one 
can then identify the driving function, and verify that the obtained maps g t 
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satisfy Loewner's equation. The form of the normalized maps and of the 
driving function already appears in appendix A of 0, but for completeness 
we present a short account of the derivation of these results below. A simple 
rescaling argument then allows us to make the connection with section 0] and 
prove that the square-root forcing we considered there generates a half-circle 
of radius \/2 when k = 9/2. 

To obtain the normalized conformal maps, we start by constructing the 
map ips of figure |H1 as the composition of the three maps depicted there. The 
explicit form of this map is 

M*) = \ a ' + f ; ( i Z \ for^G^\ 7 [0, S ] (52) 
aj a s - f s (z) 

where a s := cos(s/2) _1 and the function f s (z) is given by 



with R s : = tan(s/2). We can now expand ip s for large z, and identify the 
capacity c(s). A simple translation will then give us the map satisfying 
Expanding ip s {z) for large z up to order 1/z gives 

il> s (z) =z + {2- 2a; 2 ) + (1 - a; 4 )/z + 0(1/ z 2 ). (54) 

Therefore, the map ip s + 2a~ 2 — 2 has the expansion of equation (0), with 
capacity c(s) = 1 — a~ 4 . The proper time reparameterization that makes the 
capacity equal to 2t is given by 

2t(s) := 1 - aj 4 = 1 - cos 4 (s/2) for s E [0, tt) (55) 

which has the inverse 

s(t)/2 = arccos ((1 - 2t) 1/4 ) for t E [0, 1/2) . (56) 

The maps g t are then given by gt(z) = tp s ( t )(z) + 2a~?L — 2. More explicitly, 
this can be written in the form 



2{z - I) 2 + 4z(l - 2tfl 2 + 2{z + l) v /(^ + l) 2 -4z(l-2t) 1 /2 

(57) 
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We know from figure El that this map takes the point 7(s(t)) onto the image 
3a~^ — 2, which implies that the driving function must be 

£(t) = 3(1 -2t) 1/2 -2. (58) 

Indeed, one can check by differentiating (J57)l with respect to time, that g t 
satisfies Loewner's equation with this driving function. 

To make the connection with section HJ we note that a rescaling of the 
form (j5J) with a = 1/ \/2 turns the driving function into a function of the 
form of equation ()26|) with k = 9/2, up to a translation. We conclude that 
this particular square-root forcing produces a line of singularities that is a 
half-circle of radius y/2. 

We conclude this section with a discussion of its relation to an old paper 
by Kufarev [6.. Kufarev considers the equivalent of Loewner's equation (J2J) 
in the unit disk, and gives an explicit solution of this equation for negative 
times. A half-plane version of Kufarev's example is obtained as follows. 
Suppose that we set £(— s) = 3v^2s for s > 0. Then it can be verified that 
the solution of (J2J) at time — s is the normalized map that takes Ti, onto Ti 
minus the half-disk of radius y/2s centered at 2y2s. Observe that for s = 1/2, 
up to a translation, this solution is just the inverse of (|H7jl at time t = 1/2. 
This is not a coincidence because generally, if one sets £(— t) := £(T — t) and 
solves Loewner's equation for both and <?_t, the two results are related 
by g-T = 9t 1 - 



6 Multivalued Forcing 

In this section we look at a special case of the more general version of 
Loewner's equation (appearing for example in 



dg t _ 
dt J g t — x 



dfi t (x) (59) 



where /it is a measure on the real line that can be time- dependent. Here we 
will take ji t to be time-independent and assigning a mass pj to the values £j 
such that YlijVj — 1- Then the equation for g takes the form 

dt tfg-tj 1 ; 
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-2-10 1 2 

Figure 7: A forcing that is a superposition of the values ±1 produces two lines 
of singularities that approach each other as time goes on. This illustration 
shows the two lines up to time t = 10. 

This can be seen as a case where the forcing £ (£) is a superposition of values, 
and we also believe that this can be described as a case where the forcing 
makes rapid (random) jumps between the values 

Equations like (|6U|) are generally easily integrated. Take the specific case 
in which the possible values of £ are ±1, which are taken on with equal 
probability. Then equation ()60|) becomes 

dg_ = 2g 
dt g 2 -l 

which then integrates to give 

g 2 /2-\ng = 2t + z 2 /2-\nz. (62) 

There are two traces, starting respectively from z — ±1 and then moving 
upward towards zoo. The traces are found, as before, by setting g equal to 
the forcing. Thus the right-hand trace obeys 

z c+ {t) 2 -2\nz c+ {t) = 1-4*. (63) 

Again we can obtain an explicit expression for the trace, like we did for 
the case of linear forcing in section EJ and for the critical curve in section |3J 
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(61) 



After substituting z c (t) = r t exp(i<p t ) and splitting the previous equation in 
real and imaginary parts, we obtain 



where <ft increases with time from the value 0o — to = ir/2. Note that 
we have filled in the results for both traces, using the fact that they must 
be symmetrically placed about the imaginary axis. We can also compute the 
asymptotics. For small values of t we obtain 

z c± (t) = ±l + zv / 2tT^ + 0(t 3/2 ) (65) 
while for large values of t the traces behave like 

z c± (t) = V(4*-l)-ln(4t-l)± ^ + 0(lnt/t 3 / 2 ). (66) 

We see that as the traces move upward, they also approach one another, as 
is shown in figure [7| 
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